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Resonant tunnelling through an Anderson impurity is investigated by employing a new perturbation 
scheme at nonequilibrium. This new approach gives the correct weak and strong coupling limit in 
U by introducing adjustable parameters in the self-energy and imposing self-consistency of the 
occupation number of the impurity. We have found that the zero-temperature linear response 
conductance agrees well with that obtained from the exact sum rule. At finite temperature the 
conductance shows a nonzero minimum at the Kondo valley, as shown in recent experiments. The 
effects of an applied bias voltage on the single-particle density of states and on the differential 
conductances are discussed for Kondo and non-Kondo systems. 
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Electronic transport through an artificially fabricated 
quantum dot (QD) is of considerable current interest (for 
a review see e.g. Ref. [j] and references therein). In ad- 
dition to the Coulomb blockade and ordinary resonant 
tunnelling through a discrete electronic level,.-, a- novel 
Kondo-assisted tunnelling has been predicted!] El This 
phenomenon-has been verified in single-electron transis- 
tors (SET)Lli3 as well as in nanometer-sized metallic con- 
tacts!] In the SET structures several parameters can be 
controlled by external gates. This opens the possibility 
of studying a variety of electron correlation effects which 
are not available in the traditional Kondo problem of di- 
lute magnetic impurities in a metal host. For example, 
the level position of the quantum dot and the coupling 
to reservoirs xan be tuned in the single electron transis- 
tor structureOu, so that a tunable Kondo effect could be 
observed. Further, by applying an external voltage be- 
tween two reservoirs, it is possible to study a new kind 
of Kondo effect out of equilibrium. 

Transport through an Anderson impurity out of 
equilibrium has been studied theoretically by differ- 
ent methods including, the perturbation treatment in 
Coulomb reaulsion Un a and the non-crossing approxima- 
tion (NCA)El. The second order perturbation treatment 
gives reliable results for the symmetric case (2eq + U = 0, 
Eq being the level position of the impurity). However, it 
is known that away from the symmetric case second or- 
der treatment does not reproduce the correct low and 
high energy limitsJj For the limit of infinite U with weak 
hybridization, NCA provides a quantitatively accurate 
results for the Kondo physics except for very small tem- 
peratures. Meanwhile, NCA has some drawbacks in the 
low-energy region: it fails to satisfy Fermi-liquid rela- 
tionsEH and, as a consequence, overestimates the conduc- 
tance due to the Kondo resonances at low temperature. 
Although the qualitative physics of electron transport 
through the Anderson impurity is now well understood, a 
theory which can describe properly the wide range of pa- 
rameters, including the physics of Kondo resonant trans- 
port, charge fluctuation and the empty site limit, is still 
missing. 



In this paper, we investigate the transport through 
an Anderson impurity by means of a new perturbation 
treatment which gives the correct weak and strong cou- 
pling limits in a self-consistent manner. Therefore, one 
can deal with an interpolation scheme between two lim- 
its. This treatment has been suggested by Kajueter 
and Kotliarliil for an application to the Hubbard model 
in the high dimension limit. We extend this idea to 
study the transport through an Anderson impurity in 
a non-equilibrium regime. We have found that the zero- 
temperature linear-response conductance obtained from 
the density of states agrees-jwell with those obtained from 
Langreth's exact relation.EJ Such good agreement has not 
been obtained before by the previous theories. Starting 
from this good agreement, we discuss the effects of finite 
temperature and of finite voltage on transport. 

The resonant tunnelling through a single quantum 
state can be described by the Anderson impurity model 
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where represents the single particle energy in the reser- 
voir a(— L, R) with their chemical potential difference 
being the applied voltage, that is — Hr = e V ■ The 
parameters £o, U, and ij? denote the single particle en- 
ergy in the QD, the Coulomb repulsion, and the coupling 
between QD and reservoir states, respectively. 

In the wide-band limit of the reservoirs, the current 
formula can be written in the formll3'E3 
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where T(u>) = Tl(u})Tr(u>)/T(u)) with T a (u>) — 
nJ2k \tk\ 2 3( UJ ~ £ ka) being the coupling strength between 
the QD level and the lead a, and T(u) — r^(o;) + Tr(lj). 
f a (u) = l/( e /3(--AO + i) and p a (fij) = -±Im G CT H are 
the Fermi function of lead a and the spectral density of 
states (DOS) of the electron in the QD, respectively. 
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In calculating the Green's function 
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A a (uj) — J2kea l*fc| 2 /( w — e fc)' we start with the ansatz 
for the self-energy 



E(w) = Un + 



1 -6E( 2 )(w)' 
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where n is the occupation number of the QD level which 
should be determined self-consistently, and X^ 2 -* is the 
second order self-energy in U. Note that each Green's 
function line in the Y,( 2 \ui) diagram is given by 
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The parameter a in Eq. (Q) is determined from the con- 
dition that the self-energy has the exact behavior at hiah 
frequencies, and b is determined from the atomic limit.EJ 
Both conditions lead to the expressions 



and 
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where uq is a fictitious particle number obtained from 
Gq(uj). Hence, these two parameters give the correct 
weak and strong coupling limits in any region, including 
Kondo, charge fluctuation, and even-number site limit. 

The linear-response conductance G = dl /dV\v=o can 
be obtained from the spectral density of states according 
to the current formula Eq.(||). At zero temperature, it 
reduces to 



G 



2e 2 



47rf(0)p CT (0). 



(8) 



In Fig. [l] we show the zero-temperature linear-response 
conductance as a function of the dot level for symmetric 
coupling Tl(0) — Tfl(0). For this and all the other figures 
we have chosen U = 7.27T(0). A parabolic form of T(u>) 
centred at ui = with bandwidth W — 16.62r(0) is used 
and the energy scale in the x-axis is normalized to r(0). 
According to our results (see Fig.|l|), G < 2e 2 /h in the 
Kondo limit, which implies the Kondo-assisted transmis- 
sion. At zero temperature, the linear-response conduce, 
tance may also be deduced from th^Friedel-Langrcthtil 
sum rule, which leads to the relational 
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FIG. 1. Linear-response conductance at zero temperature 
as a function of the dot level position. The conductance is ob- 
tained from the DOS (solid line) and from the Langreth's (ex- 
act) relation (dashed line). The dot-dashed line denotes the 
occupation number of the QD. U — 7.27r(0) and a parabolic 
form of r(w) is used in this and all the other figures. 

In Fig. one can see that the conductance obtained 
from the density of states presents good agreement with 
those obtained from the sum rule. It is important to 
note that such agreement could not be achieved by other 
treatments such as thq-ordinary 2nd order perturbation 
theory or by the NCAILa 
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FIG. 2. Linear response conductance at T = (solid 
line), T = 0.5r(0) (dashed line) and T = T(0) (dot-dashed 
line). The conductance is symmetric around eo = —U/2 and 
it has minimum at this value for finite temperatures and a 
flat maximum at that point for T = 0. 

The linear-response conductance at finite temperatures 
is displayed in the Fig. @. The conductance shows a 
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strong temperature dependence in the Kondo limit, while 
it is almost temperature-independent for Eq > and 
£o + U < 0. The suppression of the conductance in the 
Kondo limit at finite temperatures is related with the re- 
duction of the density of states near the Fermi level. As a 
result, the conductance has two maxima at certain points 
in —U < sq < 0, and a minimum at eo = —U/2. Note 
that the conductance is symmetric around Eq = — U/2 
because of electron-hole symmetry of the model and has 
a minimum value at this energy for T ^ 0. To our knowl- 
edge, this is the first quantitative result which shows the 
nonzero conductance minimum due to the Kondo effect. 
Due to the Kondo resonance at finite temperature, the 
peak spacing for an odd-number QD is smaller than that 
of an even- number QD. The peak spacing between two 
conductance maxima is reduced as the temperature de- 
creases, which can be attributed to the quantum fluctu- 
ation of the number of electrons in the QD. These effec 
have been observed in recent experiments for the SET J 
The asymmetry comes from a transition between Kondo 
and non-Kondo system and is reproduced in our calcu- 
lation. The Kondo effect shows up as a nonzero conduc- 
tance between the conductance peaks for the QD with 
odd number of electrons. As the temperature increases, 
the minimum conductance at eo = —U/2 decreases and 
the positions of the maximum conductance approaches 
to the bare levels, which implies that the transport prop- 
erties are similar to the ordinary resonant tunnelling at 
high temperature. 
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FIG. 3. Zero temperature single particle DOS of the 
QD level in the presence of the applied voltage. For (a) 
eo = —U/2 and (b) eo = —0.286(7 the Kondo peaks are 
suppressed by applying the voltage, while the shape of the 
spectral DOS remains almost unchanged for non-Kondo sys- 
tems, (c) e = 0.19 U and (d) e = 0.3817. 

Let us now turn our attention to the non-equilibrium 
situation. For simplicity, we assume a symmetric volt- 
age drop, that is, fiL = —fi>R = eV/2. Fig. shows the 



spectral DOS of the QD at non-equilibrium for different 
values of eo- As one can easily find, the effects of the ap- 
plied voltage are very different for the Kondo systems ((a) 
and (b)) and for the non-Kondo systems ((c) and (d)). 
In Fig. (a) and (b) one can see that the resonance in 
the Fermi level is suppressed by applying the voltage, be- 
cause the external voltage produces inelastic scattering of 
quasi-particles located between the two different chemi- 
cal potentials. For the symmetric case (2e + U = 0), 
the spectral weight in the Fermi level is transfered to 
the satellite peaks. In the Fig. || (b), we investigate the 
asymmetric case. Here, the reduction of the Kondo peak 
is less pronounced than in the symmetric case. 

The DOS for e > 0, Fig. | (c) and (d), is insensitive 
to the external voltage. This is very different from the 
behavior of the Kondo system. In this limit, the trans- 
port properties are similar to the noninteracting system 
and the differential conductance has a weak dependence 
on the temperature and on the external voltage, as we 
discuss below. 
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FIG. 4. Zero temperature differential conductance as a 
function of the external voltage. Pronounced zero bias max- 
ima are shown in the Kondo limit (eo = —U/2 (solid line), 
—0.286(7 (dotted line)), while flat minima can be shown 
in the empty site limit (eo = 0.19(7 (dashed line), 0.38(7 
(dot-dashed line)). 

The differential conductance dl/dV as a function of 
V is displayed in Fig. || for several values of eo- The 
conductance is symmetric under bias reversal, since we 
considered a symmetric voltage drop as well as symmet- 
ric coupling for two tunnel barriers. The conductance has 
a maximum at zero bias in the Kondo system, because 
the Kondo resonance is reduced by the external voltage. 
This reduction is most pronounced in the symmetric case. 
This zero-bias nasimum has been well known from previ- 
ous calculationscfla and observed in experimentsflij The 
conductance for the symmetric case increases at larger 
voltage, which is related to the ordinary resonance of the 
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QD level with the leads. In contrast, the conductance 
shows very weak zero-bias minimum for the non-Kondo 
system (e > 0) . This can be understood from the behav- 
ior of the spectral DOS in the presence applied voltage 
in Fig. H (c) and (d). If the voltage dependence of the 
spectral DOS is neglected, the conductance at zero tem- 
perature is given by 
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According to the behaviour observed in Figs. || (c) 
and (d) as well as the Eq. ( |l0| ) one can conclude that 
the differential conductance has a weak zero-bias mini- 
mum for a non-Kondo system. The experimental results 
for ihc non-Kondo valley do not show a universal behav- 
ior.B That is, both zero-bias minima and maxima have 
been shown in the non-linear conductance. This issue 
seems to be still open, while a pronounced zero-bias min- 
imum in the mixed valence limit has been predicted in a 
previous study.ES 

In conclusion, we have studied resonant tunnelling 
through a quantum dot by means of a new perturba- 
tion treatment which correctly takes into account the 
weak and strong coupling limits. In our study, we have 
found that the zero-temperature linear-response conduc- 
tance agrees well with that obtained from the exact sum 
rule. At finite temperature the conductance has been 
shown to have a nonzero minimum at the Kondo valley 
and the minimum conductance decreases with increas- 
ing temperature. We have shown that the Kondo reso- 
nance is reduced by applying a small voltage. This leads 
to a pronounced zero bias maximum of the differential 
conductance as a function of the voltage. On the other 
hand, the density of states is insensitive to the voltage 
in a non-Kondo system, and accordingly weak zero-bias 
minima occur in the differential conductance. 
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